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7V-methyl-D-aspartate (NMDA) receptor activation involves a dynamic series of structural rearrangements initiated 
by glutamate binding to glycine-loaded receptors and culminates with the clearing of the permeation pathway, 
which allows ionic flux. Along this sequence, three rate-limiting transitions can be quantified with kinetic analyses 
of single-channel currents, even though the structural determinants of these critical steps are unknown. In inactive 
receptors, the major permeation barrier resides at the intersection of four M3 transmembrane helices, two from 
each GluNl and GluN2 subunits, at the level of the invariant SYTANLAAF sequence, known as the lurcher motif. 
Because the A7 but not A8 residues in this region display agonist-dependent accessibility to extracellular solutes, 
they were hypothesized to form the glutamate-sensitive gate. We tested this premise by examining the reaction 
mechanisms of receptors with substitutions in the lurcher motifs of GluNl or GluN2A subunits. We found that, 
consistent with their locations relative to the proposed activation gate, A8Y decreased open-state stability, whereas 
A7Y dramatically stabilized open states, primarily by preventing gate closure; the equilibrium distribution of A7Y 
receptors was strongly shifted toward active states and resulted in slower microscopic association and dissociation 
rate constants for glutamate. In addition, for both A8- and A7-substituted receptors, we noticed patterns of kinetic 
changes that were specific to GluNl or GluN2 locations. This may be a first indication that the sequence of discern- 
ible kinetic transitions during NMDA receptor activation may reflect subunit-dependent movements of M3 helices. 
Testing this hypothesis may afford insight into the activation mechanism of NMDA receptors. 



INTRODUCTION 

NMDA receptors represent one of four classes of iono- 
tropic glutamate receptors (iGluRs) responsible for fast 
excitation in central neurons. Despite overall common 
tetrameric structures, each class generates currents with 
characteristic waveforms and fulfills nonoverlapping 
physiological roles. The structural details that support 
class-specific functions are largely unknown (Traynelis 
etal., 2010). 

Like all iGluR subunits, NMDA receptor subunits 
(GluN) have modular architecture. Extracellular resi- 
dues form two genetically separable N-terminal (NTD) 
and ligand-binding (LBD) domains (Wo and Oswald, 
1995). In tetrameric receptors, eight such modules 
adopt globular shapes and stack as discernible NTD 
and LBD layers (Sobolevsky et al., 2009). Each LBD 
connects through three short linkers to the Ml, M3, 
and M4 transmembrane helices but not to the interven- 
ing M2 segment, which loops only halfway through the 
membrane from the cytoplasmic side. The permeation 
pathway is lined on its extracellular half by M3 residues 
and on its cytoplasmic half by M2 residues, which also 



form the narrowest portion of the open pore (Kuner 
et al., 1996; Beck et al., 1999; Sobolevsky et al., 2002) . 

Arguments based on structural homology and pat- 
terns of accessibility to extracellular reagents support 
the hypothesis that in iGluRs, the agonist-controlled 
gate resides at the intersection of M3 helices. In func- 
tional iGluRs, the M1-M3 segments superimpose well 
with the pore structure of the KcsA channel, whose 
ligand-dependent constriction is made by the crossing 
of four M3-like helices (Doyle et al., 1998; Jiang et al., 
2002; Armstrong, 2003; Kuner et al., 2003; Sobolevsky 
et al., 2009). In the corresponding region, iGluRs have 
a conserved nine-residue stretch, SYTANLAAF, and the 
A8 of this motif corresponds to the bundle crossing of 
KcsA (Chang and Kuo, 2008). Consistent with a loca- 
tion at or internal to the agonist-controlled gate, gluta- 
mate increases the accessibility of extracellular reagents 
to A7, whereas the more superficial A8 is accessible re- 
gardless of agonist presence (Beck et al., 1999). Func- 
tional evidence from A7- and A8-substituted receptors 
revealed stark differences in function for these two posi- 
tions in NMDA versus non-NMDA receptors. 
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When present in 8-type iGluRs, A8T substitutions cause 
gross gating modifications that result in cerebellar neu- 
rodegeneration and ataxic behaviors in the eponymous 
Lurcher mice (Zuo et al., 1997). Similarly, when intro- 
duced in AMPA- or kainate-type iGluRs, the lurcher sub- 
stitution (A8T) produces substantial constitutive activity 
and slows desensitization and deactivation kinetics of 
agonist-activated whole-cell currents (Schwarz et al., 
2001). In contrast, in NMDA receptors, lurcher muta- 
tions cause only mild kinetic changes whether present 
on GluNl (Nl) or GluN2 (N2) subunits (Zuo et al., 1997; 
Kohda et al., 2000; Hu and Zheng, 2005a,b). Instead, 
A7Y mutations (lurcher-like) and other bulky side chains 
introduced at A7 sites render receptors constitutively 
open (Jones et al., 2002; Yuan et al., 2005; Blanke and 
VanDongen, 2008). 

If indeed the adjacent A7 and A8 residues have differ- 
ent relationships to the activation gate of the NMDA re- 
ceptors, substitutions at these positions should affect 
the gating reaction with distinct kinetic mechanisms. To 
test this corollary, we determined the gating mechanism 
of NMDA receptors carrying lurcher (A8T) or lurcher- 
like (A7Y) mutations in Nl or N2A subunits. We found 
that at A8, T had no effect and Y caused receptors to 
close faster, whereas at A7, Y acted as a kinetic break on 
channel closure and shifted the entire equilibrium dis- 
tribution such that the association and dissociation rate 
constants for glutamate were slower, and a substantial 
fraction of channels was open in the absence of agonist. 
Our results are consistent with the hypothesis that A7 resi- 
dues of Nl and N2 subunits form the agonist-controlled 
gate of NMDA receptors. 

MATERIALS AND METHODS 

Cell culture and receptor expression 

RatGluNl-la (U08261) (Nl) orGluN2A (M91561) (N2) was sub- 
cloned into pcDNA3.1(+), and the following substitutions were 
introduced in Nl: A653T (N1 8T ), A653Y (N1 8Y ), or A652Y (Nl w ); 
and in N2: A651T (N2 8T ), A651Y (N2 8Y ), or A650Y (N2 W ). These 
constructs were transiently expressed in HEK293 cells (CRL-1573; 
American Type Culture Collection) along with GFP in a 1:1:1 
ratio to produce the following NMDA receptors: N1/N2 (wild 
type [WT]), N1 8T /N2, N1/N2 8T , N1 8Y /N2, N1/N2 SY , N1 7Y /N2, or 
N1/N2 W . 

Electrophysiological recordings and analyses 

Stationary single-channel currents were recorded from cell-attached 
patches with glass pipettes filled with extracellular solution con- 
taining (in mM): 150 NaCl, 2.5 KQ, 1 EDTA, 10 JV-(2-hydroxyethyl) 
piperazine-iV-(4-butanesulfomc acid) (HEPBS), and agonists (Glu, 
Gly) as indicated, adjusted to pH 8.0 (NaOH). Currents were am- 
plified and low-pass filtered (10 kHz) with an Axopatch 200B 
(Molecular Devices), digitized (20 kHz) with an A/D board (PCI- 
6229; National Instruments), and acquired and stored as digital 
files with QuB software. Files that reflected activity from a single 
active receptor throughout were selected for processing. Idealiza- 
tion and modeling were done with QuB functions as described in 
detail by Kussius et al. (2009). 



Overall kinetic parameters and gating rate constants were esti- 
mated by fitting state models to the entire sequence of events in 
each file; results are presented as the average value for each data- 
set (mean ± SEM). Statistical differences were evaluated with 
Student's t test or with one-way ANOVA and were considered sig- 
nificant for P < 0.05. Agonist binding and dissociation rate con- 
stants were estimated by globally fitting kinetic models directly 
to data obtained at several agonist concentrations, as described 
in detail by Popescu et al. (2004). Values were rounded to the 
90th percentile. 

Macroscopic responses were recorded as whole-cell currents 
with glass pipettes containing intracellular solution (in mM): 135 
CsF, 33 CsOH, 2 MgCl 2 , 1 CaCl 2 , 10 HEPES, and 11 EGTA, ad- 
justed to pH 7.4 (CsOH) . Cells were perfused with extracellular 
solution containing (in mM): 150 NaCl, 2.5 KC1, 0.5 CaCl 2 , 0.01 
EDTA, 10 HEPBS.and 0.1 Gly, adjusted to pH 8.0 (NaOH). Re- 
sponses were elicited by switching into solutions with added Glu 
(1 mM). Currents were amplified and analogue filtered (2 kHz; 
Axopatch 200B), sampled (5 kHz; Digidata 1440A), acquired with 
pClamp software (Molecular Devices), and analyzed offline. 

Simulated macroscopic currents were calculated from the spec- 
ified kinetic models, expanded as necessary with glutamate-binding 
steps. Responses were initiated by positioning all channels (ra = 
100; 10 pA/each) in the glutamate-free state C° and stepping Glu 
concentrations instantaneously from 0 to 1 mM. Responses were 
represented as time-dependent occupancies of the aggregated 
open state and were analyzed as the experimentally recorded cur- 
rents (Popescu et al., 2004). 

Energy diagrams were calculated using the rate constants speci- 
fied in each model and the relationship AG 0 = —RT(lnKeq), where 
R is the molar gas constant, T is the absolute temperature, and K<, q 
is the equilibrium constant of the transition considered, calcu- 
lated as the ratio of the forward to reverse rate constants. Barrier 
heights were calculated with the relationship E* = AG 0 + (10 — ln£ + ) 
and were arbitrarily illustrated at equal distances along the reaction 
coordinate (Kussius and Popescu, 2009) . 

Online supplemental material 

Table SI shows closed kinetic components for lurcher and 
lurcher-like N1/N2A receptors. Table S2 lists open kinetic 
components for lurcher and lurcher-like N1/N2A receptors. 
Fig. SI illustrates rate constants for NMDA receptors with 
lurcher and lurcher-like mutations estimated with 5C40 ki- 
netic models. The online supplemental material is available at 
http://www.jgp.org/cgi/content/full/jgp.201210786/DCl. 



RESULTS 

To test the hypothesis that A7 but not A8 of the SYTAN- 
LAAF motif forms the agonist-controlled gate of NMDA 
receptors, we set out to delineate the reaction mecha- 
nism of GluNl /GluN2A receptors carrying A8T or 
A7Y substitutions in GluNl (Nl) or GluN2A (N2) sub- 
units (Fig. 1, A-C) . We recorded on-cell currents from 
patches containing exactly one active receptor in the 
presence of maximally effective agonist concentrations 
(1 mM Glu and 0.1 mM Gly), low proton concentrations 
(pH 8; 10 mM HEPBS), and the absence of divalent cat- 
ions (150 mM NaCl, 2.5 mM KC1, and 1 mM EDTA) 
(Figs. 2 A and 3 A) . In these conditions, the recorded 
traces mark the stochastic flux of sodium across the 
membrane and faithfully chronicle receptor excursions 
into and out of ion-permeable conformations (Popescu 
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and Auerbach, 2003). Such traces inform directly about 
each receptor's unitary current amplitude and after 
kinetic analyses and modeling of dwell times can help 
define the receptor's activation mechanism (Kussius 
etal, 2009). 

A8 and A7 substitutions decreased NMDA receptor 
unitary current amplitudes 

In these records, we noted that regardless of the sub- 
unit in which they were introduced, A8T and A7Y 
substitutions resulted in significantly smaller unitary 
currents (Fig. 1 D and Table 1 ) . Previously, Kohda et al. 
(2000) reported that A8T produced macroscopic cur- 
rents with lower noise than WT and concluded that 
NMDA receptors carrying the lurcher mutation may 
have "at least one low-conductance open state" (Kohda 
et al., 2000) . Similarly, when channels with cysteine sub- 
stitutions at A7 of N2A subunits were modified with 
chemical reagents, the mean single-channel current am- 
plitudes decreased relative to the parent A7C mutant 
(Yuan et al., 2005). Here, we show that substitutions 
with natural residues at A7 or A8 of the lurcher motif 
also significantiy decreased NMDA receptor unitary cur- 
rent amplitude. In addition, we noted that when all four 
subunits contained A7Y substitutions, the amplitude of 
unitary currents was not further reduced relative to 
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Figure 1 . Lurcher-motif residues control NMDA receptor uni- 
tary current amplitude. (A) Sequence alignment of the exter- 
nal segment of M3 highlights the conserved lurcher motif; the 
residues tested here are in red (A7) and blue (A8). This color 
scheme is used in all subsequent panels. (B) Cartoon of iGluR 
subunit membrane topology. (C) Partial model of two diagonally 
situated subunits based on an AMPA receptor (Protein Data Bank 
accession no. 3KG2) and on KcsA (Protein Data Bank accession 
no. 1BL8) structures. (D) Unitary current amplitudes of NMDA 
receptors carrying lurcher (A8T) or lurcher-like (A7Y) substitu- 
tions in the Nl or N2A subunits. *, statistically significant differ- 
ences relative to WT (P < 0.05; one-way ANOVA) . 



the single-subunit substitutions (not depicted). These 
results highlight a role of SYTANLAAF residues in con- 
trolling the large conductance characteristic for NMDA 
receptors and advocate for a systematic investigation 
into the mechanism by which this control occurs. 

Aside from precise information about unitary current 
amplitudes, single-channel traces also illustrate how long 
receptors dwell in nonconductive (closed [C]) and con- 
ductive (open [O]) conformations, and reveal the pre- 
cise sequence in which the receptor moves between 
closed and open structures. We used statistical methods 
to extract this kinetic information, and we present re- 
sults for lurcher and lurcher-like mutations below. 

A8T substitutions had minimal effect on NMDA 
receptor gating 

We measured open probability (P Q ) , mean open times, 
and mean closed times from currents produced by WT 
(n = 5) and receptors carrying the lurcher mutation 
A8T. These parameters were not statistically different 
across the three datasets, regardless of whether the mu- 
tation was introduced in the Nl (N1 8T /N2, n = 12 and 
P > 0.05) or the N2A subunit (N1/N2 8T , n= 7 and P > 0.05) 
(Table 1) . For N1 8T /N2, we noticed a trend toward lon- 
ger closures (21 vs. 10 ms) and longer openings (17 vs. 
1 1 ms) , but these differences did not achieve statistical 
significance even for the larger size of this dataset. 

Next, we used statistical methods to determine the 
number of kinetically distinct conformations visited by 
receptors during gating and to estimate gating rate con- 
stants (Fig. 2, B and C). Single NMDA receptors gener- 
ate current traces whose open and closed intervals have 
complex distributions, each containing several compo- 
nents (Howe et al., 1988; Gibb and Colquhoun, 1991, 
1992). Previous studies concluded that in the condi- 
tions used here (high agonist concentration and absence 
of divalent cations) , interval durations have five closed 
components and up to four open components (Popescu 
and Auerbach, 2003; Kussius et al., 2009). Closures re- 
flect three kinetically distinct preopen states (C3_i) and 
two desensitized states (C4 and C5) (Banke and Traynelis, 
2003; Popescu and Auerbach, 2003; Auerbach and Zhou, 
2005; Kussius et al., 2009). Openings reflect the two obli- 
gatory components present in each mode (Oi and O9) , 
with additional multiplicity resulting from the distinct 
durations of O2 in each mode (O2L, O2M, and C>2h) 
(Jahr and Stevens, 1987; Popescu and Auerbach, 
2003, 2004; Magleby, 2004; Amico-Ruvio and Popescu, 
2010; Popescu, 2012). 

The records we obtained from N 1 8T /N2 and N 1 /N2 8T 
had interval distributions that were very similar to those 
described previously for WT receptors, an indication 
that, like WT, these receptors followed a basic kinetic 
mechanism with five closed and two open states, and 
also experienced modal gating (Fig. 2 B and Tables SI 
and S2). As described previously, the recorded activity 
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TABLE 1 

Kinetic attributes ofNMDA receptors with A 7 and A8 substitutions 



Mutation 





Amplitude 
/ 


P„ 


Mean open time 

MS 


Mean closed time 


n 


Duration 

mill 


Events 
xlO 6 


N1/N2 


10.0 ±0.6 


0.54 + 0.03 


9 + 1 


8 + 1 


5 


244 


1.5 


N1 8T /N2 


7.3 ± 0.3" 


0.48 ± 0.06 


17 + 4 


21 + 5 


12 


552 


2.8 


N1/N2 8T 


7.5 ± 0.5" 


0.56 ± 0.05 


11 ± 1 


10 + 2 


7 


266 


1.7 


Nl 8 */^ 


8.3 ± 0.8 


0.56 + 0.1 


9 + 2 


8 + 3 


6 


441 


3.1 


Nl/^ 8 " 


7.1 +0.3" 


0.11 ± 0.03" 


5 ± 0.4" 


56 ± 14" 


7 


198 


0.46 


Nl w /N2 


7.2 + 0.2" 


0.7 ± 0.03" 


65 ± 13" 


28 + 6 


5 


181 


0.23 


(-/Gly) 


5.9 + 0.4"' b 


0.24 ± 0.04"' b 


30 ± 2 a b 


112 + 25"' b 


5 


181 


0.19 


("/-) 


5.8 ± 0.2"' b 


0.06 ± 0.01"' b 


77 ± 8" 


1,400 ± 400"' b 


5 


370 


0.03 


N1/N2 W 


6.2 + 0.2" 


0.81 ± 0.06" 


110+ 10" 


29 + 13 


5 


190 


0.18 



WT 
A8T 

A8Y 
A7Y 



"Significant difference relative to WT; P < 0.05 (Student's f test). 

""Significant difference relative to Nl w /N2 with Glu/Gly; P < 0.05 (Student's nest). 



can be comprehensively described with a simplified 
5C10-state model, where each C represents a kinetically 
distinct closed state and O represents collectively all ob- 
served open states. In this simplified model, the rate 
constants for entry into and exit from this aggregate 
state report the kinetics of several microscopic transi- 
tions, including mode switching (Fig. 2 B) (Kussius et al., 
2009). Consistent with the mean parameters we estimated 



for A8T mutants, kinetic modeling results indicated 
that relative to WT receptors, gating changes were small 
for N1 8T /N2 and negligible for N1/N2 8T receptors 
(Figs. 2 C and SI). We noted that for N1 8T /N2, rate 
changes clustered around the state (with faster exit 
rates C2— >Ci and C2— >C^ and slower entry rates G^*— Ci 
and C9<— C4); in addition, entry into the principal de- 
sensitized state C5 was also faster (C3— >C$). In contrast, 
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Figure 2. Gating effects of the lurcher mu- 
tation A8T on NMDA receptors. (A) Con- 
tinuous 30-s segments were selected from 
minutes-long recordings obtained from one 
WT or one lurcher (N1 8T /N2 or N1/N2 8T ) 
receptor (on-cell; open is down). (B) Dwell- 
time distributions calculated from the records 
illustrated in A are overlaid with probability 
distribution functions calculated from 5C40 
models (thick) and individual closed or open 
kinetic components (thin). (C) Kinetic mod- 
els optimized by fits to the entire sequence 
of closed and open intervals in each record; 
for each postulated transition, rate constants 
are given as rounded means for each dataset 
(in s~'). *, significant difference relative to WT 
(P < 0.05; Student's t test). All states (C, O) 
represent fully liganded (2 Glu, 2 Gly) recep- 
tors. (D) Whole-cell currents recorded during 
5-s application of 1 mM Glu (gray) are over- 
laid with the trace simulated with the corre- 
sponding kinetic model in C (black, purple, 
or green). Traces were normalized to peak. 
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for N1/N2 the only rate change we observed was the 
resensitization rate C.3<— C5, which was faster. 

The deduced kinetic models predicted macroscopic 
responses that were comparable with experimentally 
measured whole-cell currents elicited with long (5-s) 
glutamate pulses (in the presence of glycine) (Fig. 2, 
C and D). For WT and N1/N2 8T , we observed a nearly 
perfect match between model predictions, which were 
deduced from cell-attached recordings, and experi- 
mentally measured traces, which were obtained from 
whole-cell recordings. For N1 8T /N2, the simulated re- 
sponse overlapped with only a set of the whole-cell 
traces measured. However, we noted that the single- 
channel records of N1 8T /N2 receptors also displayed 
considerable patch-to-patch variability, and this was 
consistent with the range of whole-cell traces recorded. 
Based on these results and observations, we conclude 
that overall, the A8T substitution had no observable 
effect on NMDA receptor gating kinetics when intro- 
duced in N2 subunits, and had only a modest effect 
when introduced in Nl subunits. Kohda et al. (2000) 
reported changes in macroscopic desensitization and 
deactivation for lurcher mutants. However, when trace 
amounts of contaminant Zn 2+ were removed, as is the 
case with our recordings, these changes were substan- 
tially attenuated (Hu and Zheng, 2005a). Collectively 
with these previous reports, our results support the 
conclusion that lurcher mutations have minimal effect 
on NMDA receptor gating. 



Lurcher-like mutations redistribute NMDA receptors 
into open states 

In contrast with the small gating effects we observed for 
A8T mutations, A7Y produced robust and reproducible 
kinetic changes regardless of the subunit in which it was 
introduced. We noted that the patterns of change were 
subunit dependent, and overall, the effects were larger 
when the mutation was in the N2 subunit (Fig. 3, A-C 
and Table 1 ) . Relative to WT, A7Y channels had sub- 
stantially higher P Q . The increased activity was entirely 
caused by substantially longer openings (6-fold or 
10-fold). In both preparations, closures were also longer, 
but these differences were smaller and not statistically 
significant (approximately threefold; P > 0.05) (Table 1 ) . 
Unlike the nonadditive effect observed for changes in 
unitary current amplitudes, gating changes appeared 
to add up such that open probabilities for the double 
mutant N1 7Y /N2 7Y were in excess of 0.99 (not de- 
picted). We conclude that the observed changes in 
unitary current amplitude and channel gating are not 
correlated, an indication that these effects may arise 
by separate mechanisms. 

Analyses of interval durations revealed that in all A7Y 
records, closed interval distributions had five compo- 
nents, whereas open distributions had only three com- 
ponents, lacking the open component indicative of low 
mode (t l ) (Fig. 3 B and Tables SI and S2). Relative to 
WT, the principal open components (t m and t h ) were 
five- to sevenfold longer for both Nl w /N2 and N1/N2 7Y . 



A N1 7Y /N2 




1 10 10 2 10 3 10 4 1 10 10 2 10 3 

Interval durations (ms^ 



255* ~ 1 " 280* 
N1/N2 8Y (n = 7) 



Figure 3. Gating mechanism of A7Y 
and A8Y NMDA receptors. (A) Contin- 
uous 30-s traces produced by one Nl'V 
N2 (top) or one N1/N2 receptor 
(bottom); open is down. (B) Dwell-time 
histograms for the records shown in 
A; overlaid are probability density func- 
tions calculated with a 5C30 model 
(thick line) and kinetic components 
(thin lines). (C) Reaction mechanisms 
derived from fits to the entire event se- 
quence in each file; rate constants (s ) 
are given as the rounded mean for each 
dataset. *, significant differences rela- 
tive to WT (P < 0.05; Student's t test). 
All states (C, O) represent fully liganded 
(2 Glu, 2 Gly) receptors. (D) Whole-cell 
responses to 5-s applications of 1 mM 
Glu were recorded from multiple cells 
(gray) and are superimposed with traces 
simulated with models in C (purple and 
green). (E) Occupancy plots calculated 
from the corresponding models in C. 
(F) Reaction mechanisms derived from 
records produced by NMDA receptors 
with A8Y substitutions. 
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We conclude that for both A7Y mutants, the dramatic 
increase in mean open durations occurred by two mecha- 
nisms: a substantial increase in the duration of openings 
(T t , t m , and t h ) and an absence of low-mode openings. 

To organize this information into a coherent even 
if simplified reaction mechanism, we used the same 
5C10 scheme described above for A8T mutants. The 
rate constants estimated with this approach indicated 
that although both A7Y mutations affected almost all 
rate constants, the increased P () originated primarily 
from substantially slower closing rates: approximately 
fivefold and approximately eightfold for Nl and N2 
mutants, respectively (Fig. 3 C) . Macroscopic responses 
simulated with the deduced models predicted well the 
experimentally recorded whole-cell current relaxations 
during prolonged glutamate application (Fig. 3 D). 
This close match indicated that the kinetic models de- 
duced from single-channel traces captured correctiy 
the salient changes in activation mechanism produced 
by the A7Y mutations. Specifically, the mechanisms 
postulate that A7Y mutants had dramatically more sta- 
ble open states, resulting in considerable increased 
open occupancies (125 or 155% relative to WT) at the 
expense of preopen (C3-C1) and desensitized (C4, C5) 
states (Fig. 3 E) . Thus, consistent with their hypothe- 
sized distinct physical locations relative to the activa- 
tion gate, A7Y and A8T affected gating with distinct 
mechanisms: A7 but not A8 substitutions prevented 
gate closure. 

We considered whether the observed effects reflected 
the difference in the residues' substitutions, threonine 
at A8 versus tyrosine at A7, rather than the residues' 
position along the M3 helix. To examine whether tyro- 
sine can prevent channel closure when inserted at A8, 
we introduced A8Y substitution in Nl and N2 subunits 
and determined the gating reaction mechanism for 



N1 8Y /N2 and N1/N2 8Y receptors. Results show that 
when introduced at A8, tyrosine residues had no effect 
(N1 8Y /N2, P > 0.05) or accelerated channel closure two- 
fold (N1/N2 8Y , P < 0.05) (Fig. 3 F). We conclude that 
bulky substitutions at A7 but not at A8 stabilize recep- 
tors in open conformations. These results support the 
proposed location of A7 residues at the agonist-controlled 
gate of NMDA receptors and suggest a foot-in-the-door 
mechanism for how voluminous residues prevent chan- 
nel closure. 

When we compared the C3— >C%— >C\— >0 activation se- 
quences of the six NMDA receptor mutants for which 
we developed detailed reactions mechanisms, we were 
surprised to note patterns of change that were common 
for the Nl or the N2 mutants but different between the 
two groups (Figs. 2 C and 3, C and F). For example, the 
Nl but not the N2 mutations affected the C2-C1 equilib- 
rium for all mutations examined, whereas N2 but not 
Nl mutations affected preferentially the C3— >C? transi- 
tion. The simplest interpretation of these observations 
is that exit from the kinetic states C3 and C2 is controlled 
by rate-limiting events that reflect rearrangements at 
the level of the A7 residues within N2 and Nl subunits, 
respectively. However, because the changes in closed 
components were relatively small and the inherent vari- 
ability of single-molecule measurements reduced the 
resolution with which differences could be reliably ascer- 
tained, this hypothesis is highly speculative at this time. 

Constitutive activation mechanism of N1 7Y /N2 receptors 
The experiments described above were deliberately per- 
formed in the presence of supra-maximal concentra- 
tions of glycine (0.1 mM) and glutamate (1 mM) to 
focus on transitions experienced by liganded receptors. 
The results demonstrated that A7Y substitutions stabi- 
lized receptors in open conformations. However, they 
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Figure 4. Constitutive activity of Nl /N2 receptors. 

(A) Continuous activity (30 s) selected from recordings 
obtained from one receptor with no agonists added (—/—) 
or no glutamate added and 0.1 mM Gly ( — /Gly). 

(B) Histograms of closed events observed in the correspond- 
ing records shown in A; probability distribution function 
(thick line) and kinetic components (Ei-E 6 , thin lines) cal- 
culated with a 6CIO model; the concentration-dependent 
component areas are in red. (C) Kinetic model used to 
fit the entire sequence of events in each (—/Gly) record; 
rate constants (s -1 ) are given as the rounded means for the 
entire dataset. All states (C, O) represent glycine-bound 
glutamate-free receptors. 
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did not inform about any changes that may have occur- 
red in glutamate-binding kinetics. Similar to Blanke and 
VanDongen (2008), we also observed substantial leak 
currents when measuring whole-cell activity from 
N1 7Y /N2 or N1/N2 W receptors, an indication of consti- 
tutive activity. To investigate the mechanism by which 
A7Y substitutions produced spontaneous openings, we 
set out to delineate the activation reaction of agonist- 
free N1 7Y /N2 receptors. 

We recorded one-channel activity from N1 7Y /N2 chan- 
nels in the absence of added agonists (— /— ) and with 
0.1 mM glycine in the absence of added glutamate 
(— /Gly) (Fig. 4 A). These records revealed measurable 
P„ in both conditions: 0.06 ± 0.01 (n= 5) and 0.24 ± 0.04 
(n = 5), respectively, with no statistical change in mean 
open durations relative to activity recorded from this 
mutant in high concentrations of both glycine and glu- 
tamate (P > 0.05) (Table 1 ) . Thus, relative to recordings 
obtained in high agonist concentrations, the lower ac- 
tivities observed in these two conditions were entirely 
the result of longer closed periods, caused primarily by 
the presence of an additional closed component (E 6 ), 
which was more prominent in the (~/ — ) condition 
(Fig. 4 B) . As described previously, this additional closed 
component most likely reflects dwells in nonconduct- 
ing receptor species that have at least one vacant ligand- 
binding site (Popescu etal., 2004). 

The ubiquitous incidence of glycine contamination 
makes it problematic to ascertain the presence of gly- 
cine-free NMDA receptors and prevents us from mea- 
suring gating rate constants for agonist-free receptors. 
Instead, we determined the gating reaction of glycine- 
bound glutamate-free N1 7Y /N2 channels. For this, we fit 
the records obtained in high glycine concentrations 
and no added glutamate ( — /Gly) with models that in- 
cluded receptor states corresponding to the apo- and 
mono-liganded conformations of the WT receptors (C° 
and C M ). Note however that although in this model 



receptors transition between conformations designated 
C°, C M , C, and O, at all times the N2 agonist-binding 
sites are vacant. Results show that the unitary current 
amplitude and mean open duration were similar to those 
measured for fully liganded N1 7Y /N2 receptors, and the 
much lower P 0 in the absence of glutamate was entirely 
caused by slower progress into C3, the first liganded-like 
conformation. Once the receptor has surmounted this 
energetic barrier, it gates with rates similar to those expe- 
rienced by fully liganded receptors (Fig. 4 C and Table 1 ) . 

N1 7Y /N2 receptors have higher affinity for glutamate 

Encouraged by this result, we aimed to measure the 
agonist-binding kinetics of N1 7Y /N2 receptors. At present, 
we cannot extract kinetic information about the glycine- 
binding reaction because a model for glycine biding 
to NMDA receptors is not available. Instead, we set out 
to examine the glutamate-binding reaction to glycine- 
bound receptors. For this we collected six additional 
records from on-cell one-channel patches with 0.1 mM 
glycine present and three concentrations of glutamate 
(0.5, 1, and 3 pM; n = 2 for each) . These records displayed 
the additional closed component E 6 , whose time con- 
stant was concentration dependent, becoming shorter 
as glutamate concentration was increased. To estimate 
glutamate on- and off-rates, to these six files we fit glob- 
ally a state model that, in addition to the 5C10 sequence 
inferred from measurements in high glutamate concen- 
tration, also included glutamate-binding steps (Fig. 5 A) . 
The results revealed a 2.5-fold slower glutamate on-rate 
(km, 7 x 10 6 M _1 s _1 ) and a 25-fold slower glutamate off- 
rate (& 0 ff, 2.5 s _1 ) for glycine-bound N1 7Y /N2 receptors 
relative to the rate constants reported previously for 
WT receptors (Fig. 5 A) (Popescu and Auerbach, 2004; 
Kussius and Popescu, 2010). The rate constants mea- 
sured here indicate a 10-fold increase in the calculated 
microscopic affinity for glutamate {k oS /k on , 0.36 pM vs. 
3.53 pM for WT) . Simulations with the measured rates 
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Figure 5. Nl™/N2 NMDA receptors 
have increased affinity for gluta- 
mate. (A) Reaction schemes include 
explicit glutamate-binding steps 
appended to the [5C10]-liganded 
sequences represented in Figs. 2 C 
(WT, black) and 3 C (A7Y, blue). At 
all times receptors are fully occupied 
with glycine. Rate constants for WT 
(top) are from Popescu et al. (2004); 
rate constants for N1 7V /N2 (bottom) 
were estimated from global fits to 
records obtained in high glycine 
concentrations (0.1 mM) and three 
separate low glutamate concentra- 
tions (0.5, l,and3pM; « = 2foreach). 
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predicted a 10-fold shift in glutamate potency (EC 5 o, 
0.17 pM vs. 1.9 pM for WT) (Fig. 5 B), as reported pre- 
viously by Blanke and VanDongen (2008) for receptors 
expressed in Xenopus laevis oocytes (EC50, 0.45 pM vs. 
4.05 pM for WT) . 

Importantly, the detailed reaction mechanism we de- 
veloped affords us the insight needed to separate the 
observed change in glutamate EC50 produced by the A7Y 
mutation into its constituent components, i.e., changes 
in macroscopic affinity and changes in macroscopic 
efficacy (Colquhoun, 1998) . To this end, we constructed 
a glutamate dose-response curve using a hybrid model 
that combined the rates we determined here for gating 
of fully liganded A7Y receptors and the binding rates 
reported previously for WT receptors (Fig. 5 B) . By sim- 
ulating dose-response curves with this hybrid model, 
we determined that if the mutation had changed only 
gating, the EC50 should have been 0.74 pM, which rep- 
resents only a 2.5-fold change relative to WT receptors 
(Fig. 5 B, middle trace). This result indicates that the 
fourfold change in macroscopic affinity made the major 
contribution to the 10-fold shift in glutamate EC50 pro- 
duced by the A7Y mutation. 

N1 7Y mutation energizes the resting states 
of NMDA receptors 

The kinetic models we developed enabled us also to cal- 
culate relative fluctuations in the free energies of WT 
and A7Y receptors during activation (Fig. 6 A) . To com- 
pare these energetic profiles, we made the following 
assumptions. First, based on the observation that the 
channel open durations were similar for A7Ywhen mea- 
sured in (Glu/Gly) and in ( — /Gly), we aligned the cal- 
culated energy profiles for these two conditions at the 
open states. Second, based on the observations that the 
C3-C2 transition had similar kinetics for WT and ligan- 
ded A7Y receptors, we aligned these two profiles at the 
energy level corresponding to C3. This alignment also 
reflected correcdy the more stable open states we ob- 
served for A7Y relative to WT. Keeping these assump- 
tions in mind, the calculated energy diagram postulates 
that, in addition to stabilizing the aggregated open 
state, the A7Y mutation also energized the resting states 
(C° and C M ) of glutamate-free receptors. The free en- 
ergy contributed by the mutation was comparable in 
value to the free energy provided by glutamate binding 
to WT receptors, such that glutamate-bound WT and 
glutamate-free A7Y receptors differed by only ^2 k B T. 
In addition, regardless of whether glutamate was bound 
or not, the mutation stabilized the open state by a simi- 
lar ^2-k B T amount. This analysis implies that the struc- 
tures responsible for agonist binding and channel 
gating remain highly coupled in A7Y mutants; estimates 
the energetic contributions of the A7Y mutation to the 
resting and open states of the receptor; and illustrates 
that regardless of glutamate presence, the equilibrium 



point for A7Y receptors is shifted to the right, or fur- 
ther along the activation pathway, as compared with 
WT receptors. 

DISCUSSION 

We examined single-channel activity of NMDA recep- 
tors with substitutions at positions A8 and A7 of Nl or 
N2 subunits and determined their activation mecha- 
nisms. These data and results allowed us to identify gat- 
ing changes caused by substitutions at A8 or A7, which 
may reflect discrete locations for these residues relative 
to the activation gate (external vs. internal) ; to note gating 
changes caused by substitutions in Nl or N2 subunits, 
which may reflect subunit-specific contributions to spe- 
cific steps within the gating reaction (earlier vs. later); 
and to measure separately the effects of A7Y substitu- 
tion on microscopic glutamate affinity (10-fold lower 
Kd) from those on receptor gating (1.3-fold higher P D ). 
In addition, we observed that mutations at A7 and A8 
produced changes in unitary current amplitudes that 
were independent of the effects on gating. 

We noted that A8 and A7 substitutions had opposite 
effects on open durations. Relative to WT, Y at A8 
produced a decrease (for N1/N2 8Y ) or no change (for 
N1 8Y /N2) in P„, and this change was caused in part by 
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Figure 6. Proposed NMDA receptor activation sequence. 
(A) Free energy profiles for WT (black) and A7Y (purple) were 
calculated from the kinetic models given for WT in Fig. 5 A, con- 
stitutive Nl /N2 gating (dotted line) in Fig. 4 C, and glutamate- 
activated gating (filled line) in Fig. 5 A; the resulting traces were 
aligned based on the measured differences in open-state stabili- 
ties for the three conditions. (B) A possible sequence of structural 
changes that accompany NMDA receptor activation and its corre- 
spondence to the rate-limiting steps identified in single-molecule 
traces; for clarity, desensitization steps are omitted. Only one glycine- 
bound heterodimer is represented, and NTD, Ml, and CTD 
modules are omitted. The reaction is initiated from the glycine- 
bound resting state (C°, gray) and proceeds toward active confor- 
mations represented in purple for Nl and green for N2. 
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shorter openings. In contrast, Y at A 7 produced an in- 
crease in P Q for both Nl and N2 locations, and the change 
was primarily caused by dramatically longer openings 
(Table 1 and Figs. 2 and 3). These results are consistent 
with the hypothesized locations for A8 as external and 
for A7 as internal or coincident with the glutamate- 
controlled channel gate. This functional topology was 
previously proposed by Chang and Kuo (2008) based 
on accessibility experiments. It is also supported by 
homology models based on structures of inactive KcsA 
and GluA2 channels (Doyle et al., 1998; Armstrong, 
2003; Sobolevsky et al., 2009). Still, it remains unclear 
why, despite similar structures, non-NMDA receptor A8T 
(lurcher), which is external to the gate, impacts func- 
tion as dramatically as A7Y does in NMDA receptors 
(Kohda et al., 2000; Klein and Howe, 2004). 

Second, we noted asymmetric effects of identical sub- 
stitutions when introduced in N2 or Nlsubunits. Across 
all conditions tested, the following pattern emerged: 
substitutions in N2 slowed progression for the first gat- 
ing transition, C3— >C%, whereas the same substitutions 
in Nl affected the second step, Cg— >Ci. Although the 
changes were modest, the distinct contributions ob- 
served may indicate that the M3 helices of N2 and 
Nl subunits move sequentially. Over the past decade, 
several studies have addressed the issue of subunit- 
dependent gating. So far, results show that perturba- 
tions in the ligand-binding layer and in the linkers 
connecting LBDs with transmembrane helices have 
modest or no subunit-dependent effects on gating 
(Banke and Traynelis, 2003; Kussius and Popescu, 
2009, 2010; Talukder and Wollmuth, 2011). Overall, 
these results are consistent with the view that gating is 
initiated by the agonist-induced closure of the four 
LBD clefts followed by, at least, the repositioning of 
individual LBDs within dimers, movements in the LBD- 
transmembrane domain linkers, and the separation of 
M3 transmembrane helices (Armstrong and Gouaux, 
2000; Kussius and Popescu, 2010; Borschel et al., 2011; 
Talukder and Wollmuth, 2011). Our observation that 
mutations in M3 of N2 affected the C3— >C% transition, 
and mutations in M3 of Nl affected preferentially the 
C2— >C\ transition, may reflect that all structural transi- 
tions preceding the separation of M3 helices also pre- 
cede the C3— >C% transition and thus most likely they 
are amalgamated within the long-lived state C3 (Fig. 6) . 
Clearly, to substantiate these suggestions and to estab- 
lish the trajectory of conformational change during 
NMDA receptor gating with the detail accomplished 
for peripheral excitatory transmitters, additional, more 
focused investigations are necessary (Grosman et al., 
2000; Auerbach, 2010). 

Third, we report for the first time equilibrium open 
probabilities for agonist-free A7Y receptor (P Q , 0.06) 
and for glycine-bound glutamate-free N1 7V /N2 receptor 
activity (P 0 , 0.24), and we show that in both cases, the 



gating reaction was largely similar to that of fully ligan- 
ded receptors (Figs. 3 and 4) . These results demon- 
strate that grafting a tyrosine side chain at A7 of Nl 
subunits provided sufficient energy to gate NMDA re- 
ceptors with agonist-free N2 subunits and allowed us to 
estimate that this energy was only ^2 k B T shy of the 
energy provided by glutamate binding to WT receptors 
(Fig. 6). Also, we measured microscopic glutamate as- 
sociation and dissociation rate constants for glycine- 
bound Nl w /N2 receptors and determined that both were 
substantially decreased (2.5- and 25-fold, respectively) . 
Based on these results, we estimate that Nl w -containing 
receptors have an ^10-fold greater microscopic affinity 
for glutamate and suggest that, like in WT receptors, 
conformational changes at the LBD of N2 subunits 
(responsible for changes in glutamate affinity) remain 
highly coupled with conformational changes of the M3 
transmembrane helix of Nl subunits (responsible for 
moving the gate). This coupling may explain these 
receptors' previously reported decreased sensitivity to 
competitive antagonists and to allosteric inhibitors such 
as protons, zinc, and ethanol (Low et al., 2003; Hu and 
Zheng, 2005a; Gielen et al., 2009; Xu et al., 2012). Our 
results explain the constitutive activation of N1 7Y /N2 
receptors as a shift of the entire activation equilibrium 
toward active states, caused by increased energetic bar- 
riers to channel closing and tight allosteric coupling 
between the membrane-embedded gate and the extra- 
cellular LBDs, which causes changes in affinity. 

Last, our results demonstrate that, when introduced 
at A7 of Nl or N2 subunits, tyrosine substitutions de- 
creased unitary current amplitudes by ^30% (Table 1). 
This decrease was not correlated with gating changes, 
and also unlike gating changes, it was not accentuated 
when substitutions were introduced simultaneously in 
Nl and N2. A role for lurcher-motif residues in control- 
ling NMDA receptor channel conductance has been 
suggested previously, although it has not been quanti- 
fied, and the mechanism by which it occurs is unknown 
(Kohda et al., 2000; Yuan et al., 2005) . 

In summary, probing the activation sequence of NMDA 
receptors with perturbations in the lurcher region en- 
abled us to gather additional support for the idea that 
A7 forms part of the activation gate and to show that 
introducing a tyrosine residue at this position in the Nl 
subunit has two major consequences on receptor func- 
tion. First, it energizes resting receptors to a level close 
to that of glutamate-bound WT receptors, resulting in 
substantial constitutive activity; and second, it stabilizes 
open receptors by a mechanism that prevents gate 
closure, resulting in very high open probabilities within 
bursts. Overall, the mutation shifts the reaction equilib- 
rium toward active states, decreases both the glutamate- 
binding and the glutamate dissociation rates, and 
increases glutamate affinity by an order of magnitude. In 
addition, we observed patterns of kinetic changes that 
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were subunit specific, an indication that residues in the 
M3 of Nl and N2 may affect separate rate-limiting steps 
along the activation reaction. Collectively, these results 
offer a mechanistic and quantitative platform on which to 
build toward the broader goal of resolving the temporal 
trajectory of structural changes that make the NMDA re- 
ceptor activation reaction. 
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